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Abstract. In this contribution we review our recent numerical work discussing the 
essential role of the local cluster environment in assembling massive stars. First we 
show that massive stars are formed from low mass pre-stellar cores and become mas- 
sive due to accretion. Proto-stars that benefit from this accretion are those situated at 
the centre of a cluster's potential well, which is the focal point of the contraction of the 
cluster gas. Given that most of the mass which makes up a massive star in this model 
comes from the cluster environment rather than the core, it is important to model the 
molecular cloud environment accurately. Preliminary results of a simulation which ac- 
curately treats the chemistry and time-dependent thermodynamics of a molecular cloud 
show quantitatively similar star formation to previous models, but allow a true compar- 
ison to be made between simulation and observations. This method can also be applied 
to cases with varying metallicities allowing star formati on in primordial gas to be stud- 
ied. In general, these numerical studies of clustered star formation yield IMFs which 
are compatible with the Salpeter mass function. The only possible exception to this is 
in low density unbound regions of molecular clouds which lack very low and high mass 
stars. 



1. The Assembly of Massive Stars 



1.1. Are massive stars formed from dense cores? 



Since the observations of Motte et al. (1998), it has been noted that the mass function 



of dense cores, closely resembles that of the stellar initial mass function (eg. [Johnstone 
et al.|2000t|Testi & Sargent|1998| ). This has lead many to propose that there is a direct 
link between the two distributions, with the mass of the core uniquely determining 
the mass of the final stellar system formed from it ( |Alves et al 2007 ; Simpson et al. 
2008). In this scenario massive stars would be formed from massive pre-stellar cores. 
However, the Jeans mass in molecular clouds is typically less than a solar mass, so 
it is unclear how such cores would be supported against further fragmentation. This 
could perhaps be resolved with turbulent support, as proposed in |McKee & Tan| (|2003 ). 
However, turbulence creates density enhancements which would themselves rapidly 
collapse (Do bbs et al.|2005[ ). One way of overcoming this obstacle would be to invoke 
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the effects of radiative heating from surrounding small protostars, which would raise 
the Jeans mass locally (Krumh olz & McKee|2008| ). However this model is analytic and 
does not explicitly model the cluster environment. It does not explain why the region 
of gas surrounded by the low mass stars had previously avoided collapse, despite being 
centrally condensed and at high densities. 

The connection between the core mass function and the IMF was probed in [Smith 



et al. (2009a) using a hydrodynamic SPH simulation (without feedback) of a giant 



molecular cloud. In this paper, gravitationally bound 4 p-cores' were identified using a 
clump-finding algorithm that finds structures which are part of a common gravitational 
potential minimum. T-cores' represent the bound cores which will unambiguously 
proceed to form stars, and as such are a best-case scenario when compared to observa- 
tional cores. The cores identified in observations of real clouds suffer from a variety of 
effects, such as confusion and mass uncertainties (Kainulainen et al. 2009] |Pineda et al. 



2008[|di Francesco et aL|2007] ) that make the determination of the masses of bound ob- 



jects far more difficult. Since in this simulation the final outcome of the star-formation 
process is known, i.e. the mass that is converted sink particles ( [Bate et aL]|1995 ), it is 
possible to directly trace the link between the two populations. 

The left panel of Figure [T] shows the initial cumulative mass function of p-cores at 
the point that they first become gravitationally bound. The points represent individual 
cores and are colour coded by mass. The right panel shows the mass of the stellar 
systems formed from these cores after five free-fall times with the same colour coding. 
If the stellar mass function were determined solely by the core mass function, one would 
expect to find that points of the same colour would still be grouped together. However, 
this is not the case. The points are disordered and accretion from the local environment 
has allowed many stellar systems to grow larger than their original precursor core mass. 
Even more significantly, in this simulation there are no massive cores, yet at later times 
in the simulation massive stars are formed. There are three stars in this simulation with 
masses over 15 M Q ; so if not from cores, then how is the mass assembled? 



1.2. Where does the mass come from? 



This issue was addressed in Smith et al. (2009b). In a higher resolution version of the 



previous simulation which included an approximation of thermal feedback from the 
sinks the evolution of the clumps of gas surrounding three massive sink particles was 
traced with time. The left hand panel of Figure [2] shows the column density of a large 
clump of gas which is forming a large stellar cluster. At the centre of this region a 
massive star is forming which has a mass of 29.2 M Q at the end of the simulation. The 
right hand panel of Figure [2] shows the position of gas which will be accreted by the 
massive star in the next 100,000 yr. The mass that goes into the massive star does not 
come from a well-defined 'core', but instead is distributed throughout the larger clump 
volume. Within this volume there also exist dense cores of gas (shown in blue) but the 
material accreted by the most massive star is the more diffuse inter-core material. The 
gas is gathered towards the massive protostar by the global collapse of the large clump 
as a whole. 

Figure [3] shows a cartoon picture of the how mass becomes partitioned into stars 
as the cluster evolves. On small scales, bound cores of gas collapse and convert their 
local mass into stellar systems. However, there is an additional larger-scale collapse 
as gas falls to the centre of the cluster potential. This allows cores located at the focal 
point of the large scale collapse to continue to accrete gas beyond the point at which 
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Figure 1. Left The cumulative mass function of p-cores, in which the mass of the 
p-core is recorded at the time it is first bound. Each point represents a single core 
and is colour coded by mass. Right The cumulative mass function of the resulting 
sink particles formed from the cores after they have evolved for five free-fall times. 
Each point on the graph represents the total mass of sinks formed from the cores and 
the colour coding is kept the same as in the previous plot. If the stellar mass function 
were determined solely by the core mass function, one would expect that points of 
the same colour would still be grouped together. However, this is not the case. 
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Figure 2. Left The column density of gas in a cluster forming clump, at the centre 
of which a massive star is forming. The colour scale runs from 0.05 g cm -2 to 5 
g cm -2 . Right The final fate of the gas. The green dots show the positions of gas 
which will eventually be accreted by the massive sink (red dot). Black dots show the 
position of sinks and blue dots show the location of material in cores. The gas which 
will be accreted by the massive sinks is well distributed throughout the clump, and 
generally cores within this region will not be disrupted by the massive sink. 
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Figure 3. A cartoon picture of how the gas becomes locked up into stars in a 
globally contracting cluster. 



their original core is consumed. For example, the central star of the cluster shown in 
Figure [2] has a mass of 29.2 M Q at the end of the simulation, but its precursor core 
had a mass of only 0.67 M Q . This is different from the original formulation of the 
competitive accretion theory ( |Zinnecker||1982[ |Bonnell et al.||200l] ) in which Bondi- 
Hoyle accretion was assumed, as the dynamics of this region are dominated by the gas 
potential rather than the stellar potential, and the central star is essentially stationary. 
Therefore the version of competitive accretion theory that assumes that tidal accretion 
dominates ( |Bonnell et al.|200l] ) is a better model for the growth of the massive central 
star. However, the mass function seen in these regions still follows a Salpeter like slope 
despite the fact that the tidally-dominated scenario predicts a somewhat shallower upper 
mass function. 



2. Environmental Effects 

Given that most of the mass going into the massive stars actually comes from the dif- 
fuse cluster gas rather than dense cores, it is crucial that we capture this environment 
accurately in numerical simulations. In this section we shall look at three cases of 
previously ignored environmental physics and see if they change the upper end of the 
resulting stellar mass functions. 

2.1. Chemical evolution 

Generally in numerical studies of star formation the gas is either assumed to be isother- 
mal, or is modeled by a polytropic equation of state (e.g. |Larson||2005| |Jappsen et al. 



2005| ). However, in reality, the thermal evolution of gas within a molecular cloud is 



a time-dependent function of the cloud chemistry. Recently, detailed chemistry and 
cooling modules ( [Glover & Mac Low|2007[ |Glover et al.|2010| ) have been included in 



the Gadget 2 SPH code ([Springel 2005). Processes within the thermodynamic model 



include photoelectric emission from dust grains, H2 formation heating, excitation and 
photodissociation of H2, cosmic ray ionisation, heating of dust by the interstellar radi- 
ation field, molecular and atomic line cooling, compressional and shock heating from 
dynamical processes, and time-dependent H2 and CO formation and destruction. 

This full chemical treatment has now been used to model the evolution of a 10 4 M Q 
turbulent molecular cloud irradiated by the standard interstellar radiation field (Clark 
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Figure 4. Left Maximum density along the line of sight at the onset of star forma- 
tion. Right the mean CO fraction of the gas. The maximum density and CO fraction 
are poorly correlated. 



et. al., in preparation). Figure [4] shows the maximum density along the line of sight and 
the CO fraction of the central region of the cloud at the onset of star formation. The two 
distributions are poorly correlated, as in cold dense regions the CO has been depleted 
onto dust grains. Another difference with previous models is that average temperatures 
of up to 50 K are now seen in the more diffuse gas surrounding the star formation re- 
gions. This is due to photoelectric emission from dust grains and compressional heating 
of the gas. Nonetheless, despite these differences from previous studies, the actual star 
formation is broadly similar. This work is still at a provisional stage, but at present the 
upper end of the stellar mass function produced in the simulation is consistent with a 
Salpeter slope. The major advantage of this work is that it is now possible to directly 
compare between the simulations and real data using the actual observational quantities 
(Shetty et. al., in prep). 



2.2. Gas metallicity 

A clear example of a situation in which the upper end of the mass function is thought 
to differ from the standard Salpeter slope is that of primordial star formation. The first 
stars are thought to be very massive (see e.g. |Bromm & Larson||2004[ [GTover| 2005 ; 
O'Shea & Norman||2007| ), owing to the lack of effective coolants other than H2 and 



HD in zero metalicity gas, which is expected to inhibit fragmentation. Moreover, it has 
been assumed that the first stars are single, as the envelopes in which they form show no 
sign of fragmentation up to the point at which the first star forms (e.g. |Abel et al.|2002 



Yoshida et al.||2006[ ). This would lead to a mass function resembling a delta function. 



However in more recent works some evidence of binarity has been seen ([Turk et al, 



2009[ |Stacy et aL]|2010[ Clark et. al. in prep), although this would still not produce a 



full mass function. 

However, a recent study by |Clark et aL] ( |2010| ) has shown that in the presence of 



mild subsonic turbulence, even primordial gas is highly susceptible to fragmentation. 
Figure[5]shows the central region of a collapsing Bonnor-Ebert sphere of primordial gas 
with subsonic turbulence. The sink particles formed in the haloes have masses in the 
range 0. 1 M Q to 40 M Q and the fragmentation is dependent on the local properties of 
the turbulent velocity field, and hence is predicted to vary between star-forming mini- 
haloes. In conclusion, it is still to early to make predictions about the final shape of any 
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Figure 5. Column density images showing the state of Pop. Ill clouds after they 
have converted 10 percent of their mass (100 M Q ) into sinks. The sink particles are 
denoted in the images by the white dots, and in each case we centre the image on the 
first sink particle to form in the simulation. The level of turbulence is shown on each 
plot in terms of the initial sound speed c^. 



primordial stellar mass function, but it may resemble the solar metallicity case more 
than previously expected. 



2.3. Unbound regions 

In Section 1 we found that the process of global collapse was crucial to forming massive 
stars in dense stellar clusters. This raises the question; what happens when regions of 
the cloud are unbound and not undergoing collapse? |Clark et aL| ( |2005| ) showed that star 
formation still proceeds in globally unbound regions, but that the star forming efficiency 
is lower. Further to this, |Clark et aT] ( |2008| ) and Bonnell et. al. (2010 submitted) 
examine how the mass function differs in regions which have a low stellar density and 
are not globally bound in a simulated GMC. Figure [6] shows the sink mass function in 
regions with low and high stellar densities. The mass function in the unbound regions 
differs from the global mass function. Firstly, there are no massive stars as there is 
no global collapse (although since fewer stars form in this region, this result may also 
be expected as a consequence of random sampling of an invariant IMF). Secondly, the 
mass function has fewer low mass stars, since due to the lower stellar density there are 
fewer dynamical interactions and less tidal shear from the cluster environment during 
the protostars' evolution ( [Bonnell et al.||2008| ). However, these effects will be quickly 



obscured due to stars from the nearby bound clusters being ejected into the unbound 
regions. 



3. Conclusions 

In summary, we conclude that: 



1. Massive stars are formed from the global collapse of a larger region, rather than 
from a distinct pre-stellar core. 
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Figure 6. The sink mass functions at the end of the simulation of the total (top, 
black line) sink population and for those sinks formed in red high-density regions 
(over 100 stars pc" 3 ) and blue low density regions. The low density regions are 
unbound and do not form either high or low mass stars and thus produce an unusual 
IMF. 



2. A treatment of the chemical evolution of molecular clouds produces qualitatively 
similar star formation behaviour to previous work, but facilitates a better com- 
parison with observational data. 

3. Even in primordial gas, fragmentation can produce a mass function. 

4. The expected mass function of unbound regions in a molecular cloud is deficient 
in low and high mass stars compared to the standard IMF. 
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